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AMENDMENTS TO THE CLAIMS: 

Please replace the claims with the claims provided in the listing below wherein 
status, amendments, additions and cancellations are indicated. 

1 . (Currently Amended) A method of deriving blood perfusion indices for a region of 

interest (ROI) of a subject, the method comprising the steps of: 

administering a contrast agent to the subject during a dynamic imaging scan: 
converting signal intensity data from raw unages of the scan into contrast agent 

concentration data; 

deriving parameters from the contrast agent concentration data using at least one 
transport function that accounts for delay and dispersion of the contrast agen t, wherein 
the at least one transport function includes an arterial transport function h^(t) 
represented bv a first model through a vessel leading to the ROI: and 

calculating the blood perfusion indices from the derived parameters. 

2-3. (Cancelled). 

4. (Currently Amended) A method according to claim [[3]] 1, wherein the at least one 
transport function ftirther comprising using a second mod e l to r e prosont comprises a 
tissue transport function hs(t) represented bv a second model through the ROI. 
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5. (Currently Amended) A method according to claim [[4]] i further comprising the 
step of selecting an arterial input function AIFa(t) in the vessel leading to the ROI by 
searching pixels taken of the contrast agent concentration data. 

6. (Original) A method according to claim 5 further comprising the step of measuring 
the contrast agent concentration C(t) remaining in the ROL 

7. (Currently Amended) A method according to claim [[6]] i further comprising the 
step of representing ha(t) using a ganmia-variate function (GVF) in the first model such 
that: 



ha(0 = \ 



— (r-^,)"'e-^'-''^'^» it>t,) 
A, 



0 (^ < ^i) 



where 4 =^^''^(1+^1)5 ^(^) - ^''^^ ""dx is the Gamma function, ti is 
the time taken for the contrast agent to move from the initial measurement of arterial 
input function AIFa(t) to a vessel at the entry to the ROI, u\ and ai are related to the 
mean transit time and dispersion of ha(t). 

8. (Original) A method according to claim 7 further comprising the step of estimating 
ha(t) after deriving values for parameters ti and c\ and setting =0 using the equation: 
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9. (Currently Amended) A method according to claim [[8]] 7 further comprising the 
step of determining an estimate for the arterial input function AIFt(t) of the vessel at the 
entry to the ROI using the equation: 



where ® is the convolution operator. 

10. (Original) A method according to claim 9 further comprising the step of 
determining an estimate of blood flow Ft and an estimate of the tissue IRF Rc(t) from 
the deconvolution of: 

C(0-(F,/kH)AIF,(t)(S>Re(t) 
where kH=(l-Ha)/(l-Ht) is a correction constant taking into account different values of 
arterial hematocrit Ha and tissue hematocrit Ht since the contrast agent remains in the 
extracellular fraction of blood (plasma). 

11. (Original) A method according to claim 10 further comprising the step of 
determining an estimate for the tissue transport function he(t) from the estimated Re(t) 
using the equation: 




he(0 = - 4- ^e(0 
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12 (Cxirrently Amended) A method according to claim [[1 1]] 15 further comprising the 
step of determining a rise time and a mean transit time of he(t) in order to determine 
parameters a2 and oi by assuming t2=0, where t2, ai and oi are parameters related to 
the mean transit time and dispersion of h^(t) hgj[jQ. 

13. (Currently Amended) A method according to claim [[11]] 16 further comprising 
the step of determining a peak height and a mean transit time of he(t) in order to 
determine parameters gz and ti by assuming ai-O, where t2, a2 and 02 are parameters 
relating to mean transit time and dispersion of hg{t). 

14. (Currently Amended) A method according to claim [[12]] 4 further comprising the 
step of representing a simulat e d tissue transport function hs(t) using a GVF in the 
second model such that: 



where = o-2^°^r(l + ) , t2, 02 and are parameters related to the mean transit time 
and dispersion of hs(t) through the ROL 

15. (Original) A method according to claim 14 further comprising the step of 
estimating hs(t) using the derived values for parameters a2 and 02 by setting t2=0 using 
the equation: 




^2 ^-(f-f2)f<^2 
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16. (Original) A method according to claim 14 further comprising the step of 
estimating hs(t) using the derived values for paraineters aa and t2 by setting a2=0 using 
the equation: 



1 



0 



17. (Cxirrently Amended) A method according to claim [[15]] 14 further comprising the 
step of determining a simulated tissue IRF Rs(t) using the equation: 

R,{t) = \-\lh,{r)dr 

18. (Original) A method according to claim 17 further comprising the step of 
determining a simulated contrast agent concentration Cs(t) using the equation: 

Cs(t)-^(FMAm(t)m,(t) 



19. (Original) A method according to claim 18 further comprising the step of fitting 
the simulated Cs(t) to C(t) using a least squares method according to: 

^ = I(C(/)-C,(/))^ 
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20. (Original) A method according to claim .19 further comprising the step of 

optimising the parameters Ft, ti, oi, ai, 02, a2 and t2 by minimizing S iteratively. 

21. (Original) A method according to claim 20 further comprising the step of reducing 
the number of adjustable parameters by fixing ai=0 and t2=0, or fixing ai=0 and a2=0 
leading to five adjustable parameters. 

22. (Currently Amended) A method according to claim [[20]] 8 comprising the step of 
further reducing the number of adjustable parameters by fixing a relative dispersion, 
Pi=ai/(ai+ti), of ha(t) resulting in ai dependent on tj, leading to four adjustable 
parameters. 

23. (Currently Amended) A method according to claim [[22]] 49 further comprising the 
step of calculating quantitative blood perfusion indices from the optimized parameters 
of Ft, ti,ai,ai,a2,a2andt2. 

24. (Original) A method according to claim 23 wherein the perfusion indices include 
any one or more of blood flow, blood volume, mean transit time, arterial delay time, 
arterial dispersion time or relative arterial dispersion, tissue dispersion time or relative 
tissue dispersion. 

25. (Currently Amended) A method according to claim [[24]] 49 further comprising 
the step of repeating each previous step, apart from the step of selecting the AIF, on a 
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pixel-by-pixel basis to produce quantitative maps of the perfusion indices for further 

analysis and presentation. 
26-29. (Cancelled) 

30. (Original) A method according to claim [[29]] 5. wherein the vessel is an artery, the 
method further comprising determining a venous input function VIFa(t) from a draining 
vein to estimate an AIFa(t) where a selected artery has partial voluming, the vein being 
larger than the artery. 

31. (Original) A method according to claim 30 further comprising the step of 
determining the profile of VIFa(t) from the draining vein. 

32. (Currently Amended) A method according to claim [[31]] 50 further comprising 
the step of scaling AIFa(t) to have the same first-pass bolus peak area as the VIFa(t) te 
minimiz e partial voluming eff e ct from th e AIFft (t). 

33. (Original) A method according to claim 32 wherein the first-pass bolus peak areas 
of the AIFa(t) and VIFa(t) profiles are obtained by fitting the profiles to gamma-variate 
function (GVF) profiles respectively to remove contrast recirculation effects. 

34. (Currently Amended) A method according to claim 17 further comprising the step 
of determining a simulated tissue IRF Rs(t) in the case that the contrast agent does not 
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always remain in the vascular system, such as in a tumour in the subject in order to 
determine blood perfusion indices and permeability indices using: 



where ^j(0 - 



A. 



0 



E is the extraction fraction of the tracer in the blood stream that leaks out of the vessel 
into tissue, and the tracer clearance rate constant k=E*Ft A^e is a rate constant at which 
the leaked contrast agent diffuses back into the blood stream and leaves the tissue, F^is 
the blood flow and Ve is volume fraction of the extravascular and extracellular space 
(EES), 



35. (Original) A method according to claim [[33]] 34 wherein a permeability surface 
area product PS is determined by PS = -F, ln(l - £") . 

36. (Previously Presented) Computer program means for deriving blood perfusion 
indices for a region of interest (ROI) of a subject by directing a processor to carry out 
the method steps of claim 1 apart from the step of administering a contrast agent to the 
subject during a dynamic imaging scan. 
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37. (Original) Computer program means according to claim 36 further directing the 
processor to retrieve raw image data from the dynamic imaging scan of the subject after 
a contrast agent is administered to the subject. 

38. (Currently Amended) A system of deriving blood perfiision indices for a region of 
interest (ROI) of a subject, the system comprising: 

scanning means for providing a dynamic image scan of the subject during which 
a contrast agent is administered to the subject; 

processor means linked to the scanning means for retrieving raw image data 
from the scan; 

the processor means further: 

converting signal intensity data included in the retrieved raw image data into 
contrast agent concentration data; 

deriving parameters from the contrast agent concentration data using at least one 
transport function that accounts for delay and dispersion of the contrast agen t, wherein 
the at least one transport fimction includes an arterial transport function h^(t) 
represented bv a first model through a vessel leading to the ROI: and 

calculating the blood perfusion indices from the derived parameters. 

39-40. (Cancelled) 

41. (Original) A system according to claim [[40]] 38^ wherein the at least one transport 
function further comprises a s e cond mod e l is us e d to r e pr e s e nt a tissue transport 
function hs(t) represented bv a second model through the ROI. 
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42. (Original) A system according to claim 41 wherein the processor means selects an 
arterial input function AIFa(t) in the vessel leading to the ROI by searching pixels taken 
of the contrast agent concentration data. 

43. (Original) A system according to claim 42 wherein the processor means measures 
the contrast agent concentration C(t) remaining in the ROI. 

44. (New) A method according to claim 22 further comprising measuring the arterial 
input function AIFt(t) by identifying a further artery showing a delay relative to AIFa(t) 
and thereafter fitting the estimate for the arterial input function AIFt(t) to the measured 

5 AIFt(t) in order to optimise parameters ti and a\, 

45. (New) A method according to claim 44 further comprising determining a relative 
dispersion Pi value and applying the value to all other pixels of the same subjects 
assuming a constant relative dispersion. 

10 

46. (New) A method according to claim 45 further comprising determining a constant 
relative dispersion pi value for all subjects such that the arterial transport function ha(t) 
of claim 8 is described by variable parameter ti and constant Pi. 

15 47. (New) A method according to claim 46 further comprising accounting for delay 
and dispersion by: 
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(i) deriving an initial impulse residue function Ro(t) by deconvolution of 
AIFa(t) from C(t); 

(ii) determining ti by the maximum position of Ro(t); 

(iii) deriving the arterial input function AIFt(t) of the vessel at the entry to the 
5 ROI from ha(t) of claim 8, ti and constant Pi to determine oi; and 

(iv) determining an estimate of blood flow Ft and an estimate of the tissue IRF 
Re(t) from the deconvolution of: 

C(t) = (Ft/kH)AIFt(t) ®Re(t) 

where kH=(l-Ha)/(l-Ht) is a correction constant taking into account different 
10 values of arterial hematocrit Ha and tissue hematocrit Ht since the contrast agent 

remains in the extracellular fraction of blood (plasma); and 

(v) determining perfusion indices as mean transit time MTY=^ R^{t)dt \ 
blood flow BF=Ft and blood volume BV=BF*MTT. 

48. (New) A method according to claim 31 further comprising the step of scaling 
1 5 AIFa(t) to VIFa(t) in order to minimize partial voluming effect from the AIFa(t). 

49. (New) A method according to claim 1 further comprising the step of representing 
ha(t) using a Gaussian function in the first model such that: 



ha{t) = < 



Ai 



0 (/ < /i) 
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where ti > 0, A\ = yfl^ (Ti[\ + erf (ti/yflcTi)]/ 2 and erf{t) =-= \ e" dx is the error 

function, ti is the time taken for the contrast agent to move from the initial 
measurement of arterial input function AIFa(t) to a vessel at the entry to the ROI, oi is 
related to the mean transit time and dispersion of ha(t). 

5 

50. (New) A method according to claim 1 wherein the at least one transport function 
further includes a tissue transport function hs(t) represented by a second model through 
the ROI, the method further comprising the step of representing hs(t) using a Gaussian 
function in the second model such that: 

10 



hs{t) = < 



— e^^'-'^^'^^^^' {t>ti) 



0 < ti) 



where t2>0, A2 = 4l7tGi[\^-erf{t2l42a2)y2 and erf {t)=—i=\ e'""^ dx is the 
error function, ii and aa are related to the mean transit time and dispersion of hs(t). 

51. (New) A computer readable medium storing a program for deriving blood perfusion 
indices for a region of interest (ROI) of a subject by directing a processor to carry out 
the method steps of claim 1 apart from the step of administering a contrast agent to the 
subject during a dynamic imaging scan. 
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52. (New) A computer readable medium storing a program according to claim 36, the 
program further directing the processor to retrieve raw image data from the dynamic 
imaging scan of the subject after a contrast agent is administered to the subject. 
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